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We consider a simple-model population, whose individuals react with a certain delay to temporal 
variations of their habitat. We investigate the impact of such a delayed- answer on the survival 
chances of the population, both in a periodically changing environment, and in the case of an 
abrupt change of it. It is found that for population with low degree of mutation- induced variability, 
being "slow-reacting" decreases the extinction risk face to environmental changes. On the contrary, 
for populations with high mutation amplitude, the delayed reaction reduces the survival chances. 

PACS numbers; 87.10.-|-e, 0.2.70.Lq 



I. INTRODUCTION 

In the nowadays context of global warming and habi- 
tat destruction, there is an enhanced general interest in 
the impact of environmental changes on biological popu- 
lations evolution. Despite this, little has been done from 
a theoretical point of view, in the frame of evolution- 
ary dynamics modeling, towards a systematic approach 
of the role of various elements involved in these complex 
circumstances on the population dynamics. In a recent 
paper we investigated systematically the role of the se- 
lection pressure and mutation amplitude, as well as the 
impact of the quality and quantity of the habitat changes 
on the behavior of a single-species population. 

For simplicity and in order to extract the generic fea- 
tures, we considered the case of a periodically chang- 
ing environment, as in Refs. [H, [^, The case of an 
abrupt change in the environment was also addressed. 
The mean-field level of description of the chosen model 
allowed us to put the finger, for the first time, on the 
very origin of the emerging complex behavior of this 
highly-nonlinear system, that is the delicate interplay be- 
tween the different time-scale processes. The role of the 
amplitude and period of the environmental changes on 
the critical value of the selection pressure (corresponding 
to a phase-transition "extinct-alive" of the population) 
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was clarified. However, the intrinsic stochasticity, the 
dynamically-built correlations between the individuals, 
and the role of the mutation-induced variety in popula- 
tion's evolution cannot be appropriately accounted for at 
a mean-field level. 

A more refined level of description, which is an indivi- 
dual-based one, was therefore also considered. The main 
conclusions were that the inherent fluctuations do not 
destroy the phase transition "extinct-alive" , and the mu- 
tation amplitude strongly influences the value of the 
critical selection pressure, giving rise, in particular, to 
a diversity-induced resonance phenomenon Q. The 
phase diagram in the plane of the selection and muta- 
tion parameters was discussed as a function of the en- 
vironmental variation characteristics. In particular, an 
important aspect well-known to experimental biologists, 
see e.g. was emerging naturally, namely that a small 
amount of randomness, due to mutations, is beneficial 
for population's survival in the changing environment, 
while a too large amount definitely is detrimental to it. 
The differences between a smooth variation of the envi- 
ronment and an abrupt, catastrophic change were also 
clarified, pointing to the beneficial role of the mutation 
in ensuring species survival after a catastrophe. 

In this short paper we shall address another aspect of 
this survival problem, namely the role of the delay in 
the "reactions" of the individuals. The lagged response 
to environmental changes is a phenomenon widespread 
in nature [1, d, [13, [O, [H, However, an extensive 
theoretical analysis of its impact on population dynam- 
ics is still lacking. The role and effects of time-delay in 
biologically systems has been addressed previously in the 
context of Lotka-Volterra type of dynamics of interacting 
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species [ij], where the "delay" was included at the level 
of the coupling between the species. Here we are consid- 
ering a different problem, namely the delayed-response of 
the individuals of a single-species population to a chang- 
ing environment. Using a simple model, we shall try to 
clarify the degree and limits of validity of the commonly- 
spred belief that "a population of fast-reacting individ- 
uals has better survival chances face to changes in their 
environment" . 



II. MODEL 

We consider the same type of model as in Ref. 
namely a population of hermaphrodite individuals (i.e., 
which, although bisexual, need mating for reproduction), 
living on a two-dimensional square lattice of size L x L. 
We assume that the individuals cannot cross the borders 
of the lattice. Moreover, the lattice has a finite carry- 
ing capacity, which comes from the exclusion assumption 
that there is at most one individual in each lattice node. 

The dynamics of the population takes place at dis- 
crete time-steps and is the result of: natural selection 
(interaction with the environment), individual motion, 
mating and reproduction, as described below. 

A. Natural selection. Individual trait, time- dependent op- 
timum, fitness, delayed-response, selection pressure, ex- 
tinction probability. 

Each individual i is characterized by its trait or pheno- 
type, which for simplicity is represented here through a 
real number Zi G [0,1]. The trait is fixed once and for all 
at the birth of the individual. 

The population lives in an environment whose influ- 
ence on the individuals is encoded in the value of the 
so-called optimum, ip G [0,1], which we suppose to be 
homogeneous in space, but periodically variable in time 
Lp = pit). Moreover, we consider here the simplest pos- 
sibility, 

pit) = 0.5 + A sin (^-K ^—^f-^ e(t - iinit) . (1) 

Here A denotes the amplitude of the environmental os- 
cillation, with < A ^ 0.5, To is its period, and tinit is 
the moment of onset of the optimum perturbation; O is 
the Heaviside step function. 

The case of an abrupt change in the environment, for 
which the optimum jumps at i = ijnit from p — 0.5 to 
p = Q.h -\- A was also considered. 

An individual i "reacts" with a certain specific delay 
to the changes in the environment. This means that its 
instantaneous fitness (or "adequacy to the environment" , 
see below) at time t is determined by the value of the 
optimum at a previous time [t — Ti), 



following expression: 



Pi{t) = 1 - exp 



f,{t)^l-\z,-p{t-n)\. 
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The fitness determines the instantaneous individual ex- 
tinction probability per time step Pi{t), according to the 



where 5 is a parameter which models the selection pres- 
sure of the environment and constitutes a main control 
parameter of the system. During its life-time, and indi- 
vidual oscillates cyclically from being perfectly-adapted, 
when Zi = pit — Ti), i.e., from a minimum possible ex- 
tinction rate Pi{t) = cxp(— 5), to a worse adaptation, 
which corresponds to Zi ^ p[t — Ti) and to a larger in- 
stantaneous extinction probability, and finally to a total 
lack of adaptation, when Pi{t) = 1. The pool of adapted 
individuals changes thus at each time step. 

The choice (jS]) we made of the extinction probability 
and the implicit definition of the selection pressure pa- 
rameter S are frequently encountered in the biological 
literature, see e.g. [3]. Other choices and thus other 
ways of measuring the "selection pressure" are of course 
possible. However, most of them can be mapped one 
onto the other and/or account for equivalent qualitative 
aspects of the interaction between the individuals and 
their environment. 

The individual delayed-response time Ti is fixed once 
and for all at one individual's birth. We consider here 
a simple case when the t^'s are random variables drawn 
from an uniform distribution within an interval [0, T^]. 
The upper limit of this interval Td represents another 
control parameter of the model. It is obvious that for 
a periodic variation of the environment only the values 
Td < 2To are relevant. 

Note also that an equal delay-time for all individuals 
amounts simply to a change in the time-origin. As 
such, a mean-field level description of the population 
dynamics (which is already known as inappropriate for 
describing mutation, see [l|) will not be able to account 
for the effects of the individual delay-times on the 
global evolution of the population. We shall therefore 
focus exclusively on the individual-based numerical 
simulations. 

B. Individual motion. 

An individual can move to its surroundings, and the 
simplest possibility that we shall adopt hereafter is a 
random-walk. Namely, in one time step the individ- 
ual jumps on the lattice, from its initial location to 
a randomly chosen nearest-neighbor one (i.e., a site 
within the von Neumann neighborhood of the initial 
node), provided that the chosen site is empty, and that 
it lies within the boundaries of the system. If none 
of the four first-neighbor nodes is empty, then the in- 
dividual cannot move, and thus cannot mate (see below). 

C. Mating and reproduction. Heredity and mutation. 
Suppose an individual i reaches a destination node. If 
there are other individuals ( "neighbors" ) in the nearest- 
neighborhood of this destination site, then the individual 
"z" choses at random one of these neighbors, call it "j". 



3 



for mating. The pair of individuals i and j may then 
give birth to offsprings, which are placed at random on 
the empty nodes of the joint nearest-neighborhoods of 
the two parents (that counts 6 sites); therefore, the max- 
imum number of offsprings of one pair of parents equals 
6. If there is no room in this neighborhood for putting 
an offspring, then this one is not born. 

The trait of a progeny k coming from parents i and j is 
determined by the parents' traits (heredity), but it can 
also present some "variations" due to different random 
factors, such as recombination, mutations, etc. We shall 
assume that 

Zk ^ ■^{zi + Zj) + nik , (4) 

where nik represents these variations. It brings diversi- 
fication into the phenotypic pool of the population and 
we call it conventionally mutation. For simplicity, we 
shall admit that ruk is a random number, uniformly dis- 
tributed in the interval [—M, M], where < < 1 is 
called hereafter the mutation amplitude and is a control 
parameter of the system [l^. Moreover, if Eq. ^ leads 
to Zfc > 1 or Zfe < 0, then one "renormalizes" it by re- 
setting Zk to 1, respectively 0, which means simply that 
the trait of the individuals cannot overcome some fixed 
limits. This choice (|4|) for the trait of an offspring is often 
made in the biological literature [15]. 

The population dynamics is thus driven by two main 
"forces" that are acting, to some extent, in opposite 
directions: selection and mutation, characterized, re- 
spectively, through the values of the control parameters 
S and ^A. Selection, combined with heredity, tries to 
bring the average trait close to the optimum, while mu- 
tation introduces diversity in the individual traits, and 
thus is broadening the distribution of the population's 
traits. 

The Monte-Carlo simulation algorithm considers the 
individuals distributed on the lattice nodes, the initial 
condition being represented by their positions, the pre- 
scribed values of the individual traits and delay-times. 
The initial N{0) individuals are randomly-distributed 
with a mean concentration c(0) = iV(0)/L^, and their 
individual traits are randomly assigned from an uniform 
distribution between and 1. 

The individuals are evolving, at discrete Monte-Carlo 
time steps (MCS, defined hereafter), according to the 
stages ^"Cof the dynamics as described above, namely: 
^. At a given time t an individual i is picked at random, 
and its extinction probability Pi{t), corresponding to 
one MCS, is determined according to Eqs. ([2]) and ([3]). 
Then a random number r is extracted from an uniform 
distribution in the range [0, 1]; if r < Pi, the individual 
dies, otherwise it survives. 

B. If it survives, the individual i jumps at random to 
one of the empty nearest-neighbor nodes on the lattice. 

C. Then it possibly mates and produces offsprings. 



If at the time t there are N{t) individuals in the system, 
then the above steps A-C are repeated N{t) times; this 
constitutes one MCS, the unit-time of the simulations. 
Afterwards, the time is advanced by one step, t t + 1, 
and the above algorithm is repeated. 

As a last remark on the model, it is known on general 
backgrounds [Tg^ that the system size is playing a certain 
role on the location of the phase transitin point, as well 
as on its "sharpness". We used for all our Monte Carlo 
simulations a system of 100 x 100 lattice sites, for which 
we had shown previously, see Ref . [ij , that the qualitative 
features of the phase diagram are practically not affected 
by finite-size effects. 



III. RESULTS 

For a periodic oscillation of the optimum we investigate 
the temporal evolution of a population starting from a 
given initial concentration c(0). Depending on the char- 
acteristic parameters, the population can evolve, on the 
average, either to an "alive phase", for which its con- 
centration is actually oscillating periodically, with period 
To/2, around a nonzero mean value, or can get extinct 
after a transient period of time. In our previous paper [l[ 
we investigated in detail the phase diagram "extinct- 
alive" of the population in the plane of the control param- 
eters S and A4 , for different values of the characteristics 
A and Tq of the optimum oscillation. The same type of 
phase diagram was also constructed for the case of an 
abrupt jump of the the optimum. 

Our principal concern in this paper is to determine 
how the delay in the individual response to the changing 
environment - i.e., the value of the control parameter 
- affects the phase diagram extinct-alive of the system. 
We performed extensive simulations for various range of 
parameters and the main results are illustrated in Fig. [T] 

One notices several interesting features exhibited by 
these figures: 

a) Consider first the "intermediate" values of Tq for 
which, as described in Ref. [H for the no-delay case, one 
encounters the diversity-induced resonance phenomenon, 
i.e., the "peak" in the phase-diagram illustrated in the 
upper and middle panels of Fig. [1] Then: 

(i) For small values of the mutation amplitude Ai, the 
existence of a delay in the response of the individuals 
to environmental changes (i.e., Td ^ 0) is increasing the 
survival chances of the population. The diversity related 
to the randomness in the response of the individuals 
can contribute to the appearance of a larger pool of 
well-adapted individuals and is thus formally equivalent 
to an increase in the "effective" mutation amplitude, 
which is beneficial for the survival [l|. 

(ii) For large values of A4, however, adding the ran- 
domness of the delayed-response to the mutation-related 
one is leading to an even higher "effective" mutation 
amplitude. As such, the extinction risk of the population 
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FIG. 1: Phase diagram extinct (above the curve) - alive (be- 
low the curve) in the plane of the selection pressure S and 
mutation amplitude M, for different values of the delay time 
Td and of the optimum oscillation period Tq. From the up- 
per to the lower panel, To = 5000, 1000, and 50, respectively; 
the values of the other parameters are L = 100, c(0) — 0.7, 
A = 0.3, and tinit ~ 1000. The average was taken over 10 
realisations of the stochastic dynamics and the estimated er- 
rors in the value of the critical selection pressure are less than 
±0.002. 



is increased: as seen in the plots, the phase diagram 
for the populations with delayed-response (T^ ^ 0) hes 
always below the one of the instantaneously-reacting 
population {Td — 0). 

(iii) Finally, the peak related to the mutation-induced 
diversity is generally still present for the systems with 
time delay. However, in this case the randomness in 
the delayed-response can turn a part of the pool of 
well-adapted individuals into less-adapted ones, and 
thus the height of the peak is reduced as compared to 
the case of an instantaneously-adapting population. For 
large delays (like Td — 0.75 Tq in the figures) this peak 
can be even suppressed. 

b) One concludes therefore that the role of the 
delay-induced diversity is an increase in the "effective" 
mutation amplitude. As such, it can be easily predicted 
that for small values of Tq (rapid oscillations of the 
environment) the dynamics of the system will be only 
slightly affected by the delay, since it is already only 
slightly sensitive to changes in Ai. This is illustrated 
in the lower panel of Fig. [T] No diversity-induced 
peak, i.e., no optimal "effective" mutation amplitude 
is encountered in these cases, any mutation and any 
delay in response being harmful for the surviving of the 
population. 

A way to get a better insight into the reasons for this 
behavior is the monitoring of the temporal evolution of 
the pool of fittest individuals (i.e., the individuals with 
fi{t) — 1). Figure [2] iUustrates this point, for a fixed 
value of Tq and three values of ^ 0, coresponding to 
the middle panel of Fig.[TJ The upper panel of Fig.[2]per- 
tains to the region of small mutation amplitudes in the 
phase diagram, for which a delayed-response enhances 
the survival chances. The lower panel refers to the re- 
gion of the peak in the phase diagram, for which delay 
increases the extinction risk. 

One can see that for the surviving populations the 
number n(/ = 1) of the instantaneously fittest individu- 
als is oscillating periodically in time (but never reaching 
zero), while it decays (with oscillations) to zero for the 
populations that will get extinct. The pool of the fittest 
individuals is enhanced by the delay-induced diversity in 
systems with small mutation amplitude (upper panel of 
Fig. ^ and, on the contrary, it is depleted by the delayed- 
response in populations with intermediate and large mu- 
tation amplitude (lower panel of Fig. [2]) . 

Finally, we addressed also the effects of a delayed an- 
swer in the case of a catastrophic, abrupt change in the 
environment. As illustrated in Fig. [31 one encounters 
the same type of phenomena as in the case of a smooth 
variation of the optimum, namely the fact that for small 
mutation rate the largest the delay parameter Td, the 
bigger the survival chances of the population. 

In order to understand the mechanism underlying this 
behavior of the populations with small mutation ampli- 
tude A4, it is useful to follow the temporal evolution of 
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FIG. 2: Time evolution of the number n{f — 1) of the fittest 
individuals of a population for different values of the delay 
time Td. Upper panel: M = 0.025 (small mutation ampli- 
tude), lower panel: M = 0.1 (intermediate mutation am- 
plitude). The legend in the upper panel also applies to the 
lower panel.The values of the other parameters are L = 100, 
c(0) = 0.7, A = 0.3, To = 1000, and Unit = 1000, correspond- 
ing to the middle panel of Fig. [T] 



the fitness histogram "number of individuals n{f) versus 
fitness /" . This is done in Fig. |4] for two populations 
that differ only through the value of the delay parameter 
Td, such that one of them gets extinct, while the other 
one survives after the catastrophe. Before the catastro- 
phe, the histogram had an important peak at / = 1, 
and a tail (due to mutations) to low-fittnesses. After 
the catastrophe, a new peak of low-fitted individuals ap- 
pears, such that the histogram becomes bimodal. One 
notices that the existence of a larger delay time ensures 
the persistence of a sufficiently large pool of high-fitted 
individuals even after the carastrophe, and this pool will 
ensure the survival of the species till the new-born in- 
dividuals get adapted slowly, through small mutations, 
to the new environment. For a surviving population the 
histogram becomes peaked again, in the long run, around 
/ = 1. A shorter delay time T^, however, cannot en- 



FIG. 3: The phase diagram extinct (above the curve) - alive 
(below the curve) in the plane of the selection pressure S and 
mutation amplitude A4, for an abrupt jump in the value of 
the optimum, from if — 0.5 to ip = 0.8, for different values 
of the delay time Td- The values of the other parameters are 
L = 100, c(0) = 0.7, tinit = 1000. The average was taken over 
10 realisations, and the estimated errors in the value of the 
critical selection pressure are less than ±0.002. 



sure this persistence of the high-fitted individuals pool 
for a long enough time, and the population dies, since 
the adaptation through mutations is not rapid enough. 

As seen in Fig. [3l on the contrary, for large mutation 
amplitudes the larger the delay Td, the higher the extinc- 
tion risk, since, as in the case of a periodically-varying 
environment, in this case the delay-induced stochastic- 
ity adds up to the mutation, leading to an even higher 
effective mutation amplitude, which is harmful for the 
system. 



IV. CONCLUSIONS 

We considered a simple model of single-species popu- 
lation dynamics in a changing environment and we inves- 
tigated the role of a delayed answer of the individuals to 
these habitat changes. In the case of a smooth variation 
of the environment, it was found that, in general, for 
populations with small mutation amplitudes it is more 
beneficial, in terms of the survival chance, to be slow- 
reacting than to answer instantaneously to the variations 
of the environment. However, for intermediate and large 
mutation amplitudes, faster reactions are preferable to 
slower ones. In case of a very-rapidly oscillating envi- 
ronment, the rapidity of reaction influences only slightly 
the survival chances. The same type of statements holds 
true for the case of a catastrophic, abrupt jump in the 
optimum. 

As such, one has to be rather cautious with "common- 
sense" statements of the kind "a population of fast- 
reacting individuals has better survival chances face to 
changes in their environment" . Of course, more complex 
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and realistic models than the one we presented here 
are needed in order to make more detailed quantitative 
statements and reliable predictions for real biological 
systems, and to investigate further aspects of the 
intricate problem of a population evolving in a changing 
environment. 



0.7 



0.9 



200 







T, = 1000 

d 


150 




. . . T=50 

d 




t= 1100 



? 100 



50^ 




0.7 



0.8 
/ 



200 r 



150- 



? 100 - 



50- 





1000 


■ ■ ■ 


50 



t = 2100 




0.7 



0.8 
/ 



0.9 



FIG. 4: Fitness histograms: number of individuals n(/) ver- 
sus fitness / for two populations with low mutation amplitude 
M — 0.05 and different delay-response parameters Td = 1000 
(continuous line) and Td — 50 (dotted line) in case of a catas- 
trophic event. The optimum jumps at tinit ~ 1000, from 
95 = 0.5 to </5 = 0.8. The upper panel corresponds to time 
t = fOOO, just before the optimum jump. The middle panel 
corresponds to t — flOO, and the lower panel to t — 2100 
(when only the population with Td — 1000 survived). The 
other parameters are L — 100 and c(0) — 0.7. Single runs 
were considered. 
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